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ABSTRACT 

Three-dimensional trusses are designed for minimum 
weight, subject to constraints on: member stresses, Euler 

buckling, joint displacements and system natural frequencies. 
Multiple static load conditions are considered. 

The finite element displacement method of analysis is 
used and eigenvalues are calculated using the subspace 
iteration technique. All gradient information is calculated 
analytically. 

The design problem is cast as a multi-level numerical 
optimization problem. The joint coordinates are treated as 
system variables. For each proposed configuration, the 
member sizes are updated as a sub-optimization problem. 

This sub-problem is efficiently solved using approximation 
concepts in the reciprocal variable space. The multi-level 
approach is shown to be an effective technique which con- 
veniently takes advantage of the most efficient methods 
available for the member sizing problem. 

Examples are presented to demonstrate the method. The 
optimum configuration is shown to be strongly dependent on 
the constraints which are imposed on the design. 
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I. 



INTRODUCTION 



The process of optimization of structures has undergone 
important changes since its development in the early 1960s. 
Minimum weight of elastic truss structures subject to multi- 
ple loadings has been an active area of research. Attention 
has been focused on the problem of least weight, when the 
overall layout was known in advance, and when the cross- 
sectional area was the design variable. Some attention has 
been directed toward the optimum configuration of the 
structure. Design improvements in this area often exceed 
those in fixed-geometry and so shape optimization is of 
major interest. 

Pioneer work in shape optimization was conducted in 1964 
[Ref. 1] by Dorn, Gomory and Greenberg. The optimal con- 
nectivity of nodes for truss members, subject to a single 
load condition, was found and minimum weight designs were 
achieved. In their work only planar trusses were tested, 
and the process was presented as a plastic design problem 
using linear programming. 

Their work was followed by Dobbs and Felton [Ref. 2] 
who in 1969 investigated the effect of multiple load condi- 
tions on the optimum configuration of trusses through the 
use of non-linear programming methods. Again, only planar 
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trusses were considered, subject to failure by stress and 
elastic buckling. 

Later, Pedersen [Ref. 3] considered the positions of 
the joints as continuous design variables in addition to 
the areas of the bars. Stress, displacement, and buckling 
constraints were considered. Pedersen's work is significant 
because the optimization process is carried out by consider- 
ing two separate design spaces. The optimization is 
achieved by successive iterations using a gradient method 
with move limits. 

The optimization process was advanced by Vanderplaats 
and Moses [Ref. 4] who divided the design space into two sub- 
spaces, separating the area variables and the joint position 
variables. Multiple load conditions and constraints on 
stresses and Euler buckling were considered. The optimiza- 
tion was carried out alternatively between the two spaces 
until convergence was achieved. Three-dimensional indeter- 
minate trusses were designed subject to multiple loading 
conditions. This work was extended by Vanderplaats [Ref. 5] 
to include displacement constraints. 

In other research, Spillers [Ref. 6] considered statically 
indeterminate trusses. The optimization followed an itera- 
tive design where the member sizes and node locations were 
the design variables. 
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Recently, Imai [Ref. 7] treated the sizing and con- 
figuration variables simultaneously for either determinate 
or indeterminate trusses. The optimization was achieved 
using the Augmented Lagrange multiplier method. 

The design problem considered in this study is the opti- 
mum configuration of three-dimensional indeterminate 
trusses, for multiple prescribed static load conditions. 

The objective is to minimize the weight of the structure, 
design variables are the node coordinates and the member 
sizes. Constraints include stress, Euler buckling, dis- 
placement, and the natural frequency of the system. 

Some approximation concepts are introduced in order to 
reduce computational effort and to reduce the nonlinear 
characteristic of some constraints. Among these, a first 
order Taylor Series expansion is applied to approximate the 
constraints . 

The optimization proceeds iteratively in two design 
spaces: the member sizing space, where the structure is 

optimized for a fix layout, and the coordinate space, where 
the geometry is allowed to vary. In both optimization 
processes the minimum weight of the structure was maintained, 
subject to the requirements that the constraints remain 
satisfied. 

The mathematical formulation is presented in Chapter II. 
The objective function and constraints are defined in terms 
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of the design variables. The analytic gradients are also 
formulated. 

The optimization technique is discussed in Chapter III 
Several examples are presented and explained in Chapter IV 
Conclusions and recommendations for future work are 
presented. 
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II. MATHEMATICAL FORMULATION 



A. INTRODUCTION 

Several features are desirable when finite element 
methods of structural analysis are used in optimization. 
First, the number of analyses for the structure should be 
kept as low as possible. Second, the amount of gradient in- 
formation required during the design process should be 
reduced as much as possible. 

B. ANALYSIS 

1 . Static Analysis 

The initial layout of the truss, the member sizing 
and material properties (which may be different for each 
member) , a set of external loads, and support conditions 
are initially specified. 

The analysis for the stresses and deflections must 
be carried out satisfying the conditions of equilibrium of 
forces at the nodes and compatibility of deformation. If 
the material of the structure behaves in a linear manner, 
Hooke's law will establish the force-deflection relationship. 

For a truss, it is also necessary to establish the 
following assumptions before selecting the method of comput- 
ing the internal forces. The discrete element is treated 
as pin-connected, and loads and reactions are supported at 
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the joints. In this study, the weight of the members is 
not included as loads. 

The Displacement (Stiffness) method [Ref. 8] 
considers the joint displacement components as the unknowns, 
and is written in the most general form using matrix notation, 
Ku = F (Eq. 1) 

where K is the global stiffness matrix, F is the vector or 
vectors of applied loads, and u is the vector or vectors of 
displacements. Equation 1 is the set of equilibrium equa- 
tions, and is formulated such that the compatibility is auto- 
matically satisfied. 

The generality of the method is important if either 
statically determinate or indeterminate trusses are analyzed. 
The global stiffness matrix is symmetric and sparce. These 
features are used to write the code for a computer solution, 
and the matrix K is stored in compact form for efficient 
numerical solution. 

Once the displacements at every node are known, the 
internal forces or stresses are calculated by applying force- 
deflection relations. This is defined as: 





r \ 


E. D T 


u m - u kj 


- V 


v *: - v kj 




■ w kj, 



where i and j are element and load condition numbers respec- 
tively, and k and Z are node numbers associated with the 
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element i. is the stress, E^ is the Young's modulus, 

and is the length. Matrix D = f(a,y,v) contains the 
direction cosines. For brevity, hereafter the second sub- 
script is omitted and it is assumed that the stress or 
displacement corresponds to the appropriate loading condition. 

2 . Dynamic Analysis 

When system natural frequency constraints are imposed 
in the design process, the corresponding dynamic analysis of 
the structure has to be carried out. This requires the 
solution of an eigen-problem to determine the natural fre- 
quencies and normal modes. For linear elastic structures, 
the finite element approach leads to the well-known equation 
of motion, considering free vibration conditions, 

yu + gu = 0 (Eq. 3) 

where M is the global mass matrix, and u is the linear accel- 
eration. Assuming a solution of the form 

u = (Eq. 4) 

where u> is the angular natural frequency of vibration of 
the structure, and $ is the corresponding eigenvector. 

After substitution into Eq. 3, the generalized eigenvalue 
problem becomes, 

K<t> = (jJ 2 M<$ ( e< 3 • 5) 

Written in matrix form for several eigenvalues, Eq. 5 
becomes , 

m = f2 2 M3> (Eq. 6) 
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where $ is the modal matrix, and Q 2 = diag (w 2 ) is the spectral 
matrix . 

From the static analysis, the global stiffness 
ma trix is already calculated, then only the global mass 
matrix evaluation is needed. Both generalized and lumped 
mass options are coded. These are defined respectively as: 

(Eg. 7) 



ME 

!g = X,. 
1=1 



0 . A. L . 

p l l l 



21 . 



and 



ME 



?>L = X. £ 



. A.L. 

ill 



I 

. J 6x6 



(Eg. 8) 



i=l 6 

Where I is the identity matrix, the sub-space iteration 
method of Bathe and Wilson [Ref. 9] is used to solve for a 
specified number of lowest eigenvalues and the associated 
eigenvectors. The method is economically efficient for 
large problems. The mass matrix may be diagonal or banded. 
The method is well suited for re-analysis when small 
changes are made in the design. 



C. ANALYTIC GRADIENTS OF THE CONSTRAINTS 

The necessity to compute gradients of the relevent 
functions in a design optimization process arise from the 
fact that efficient mathematical programming algorithms 
reguire information on derivatives. Furthermore, approxi- 
mate methods based on a Taylor Series expansion of functions, 
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requires determination of derivatives. Based on the static 
and dynamic analyses, the gradients of forces, displacements 
and frequencies with respect to the reciprocal of the cross- 
sectional areas and the coordinate variables, are formulated 

1 . Gradient of Member Stresses with Respect to the 
Reciprocal of Area Variables 



Stress in a member is defined as 
F. 



(Eq. 9) 



also 



o . = 



K . U • 

l 



(Eq. 10) 



The partial derivative is then? 



3c. 

3X 



3X 



A. ~l 

l 



K. . 

3 

«i + a7 — 

i ~e 



UA 



(Eq. 11) 



-e ~e 

The first term of the right hand side is zero, then 



= 5. 

9X7 



4 4 ( 'i’ 



(Eq. 12) 



This can be written in. explicit form as 



3CK 

3X 0 



T 

E.D. 

Xal 

L? 



2 

3X, 



<»i» 



(Eq. 13) 



Gradient of Nodal Displacements with Respect to the 
Reciprocal of Area Variable 



Consider the equation 
Ku = F 



(Eq. 14) 
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The derivative of [K u] with respect to some variable X, 
is 




(Eq. 15) 



In this case = 1/A^ and the loads are constant and in- 
dependent of the areas (weight of the truss elements are 
ignored) , then 




(K) (u) + K 




(Eq. 16) 



or 



K 








(Eq. 17) 



Finally, then 





(Eq. 18) 



where K ^ is the inverse of K. It is necessary to compute 
the partial derivative of the global stiffness matrix K 
defined as: 



NE E . A • 

K = Z [Dj (Eq. 19) 

i=l i 

where is the matrix of direction cosines defined earlier. 
Therefore, the partial derivative with respect to the 
reciprocal of the area variables is defined as: 




NE E.Af 

X. "TV tBil 

1=1 1 



(Eq. 20) 
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where 5^ is the kronecker delta defined by 6. = 1 if i = £ 
and i = 0 if i ^ i. In practice the full matrix, ^r-(K) 

IX, a a. 

3 ' ^ 

is not actually stored. Instead the product, (K)u is 
created directly. Details of efficient gradient computations 
are given in [Ref. 10] . 

3 . Gradient of Frequencies with Respect to the Reciprocal 
of Area Variable 



Consider the eigen-problem defined by, 

[K - w?M] = 0 (Eq. 21) 

Taking the derivative with respect to the variable X [Ref. 11] 
gives 

4“ {[K - U3?M] *.} =0 (Eq. 22) 

-i 



then 



dr { '5 - “&] > + 16 - “(W 4 



<j> . = 0 (Eq. 23) 






T 

Pre-multiplying through by ^ and applying the condition of 
symmetry of the matrix gives 



il [5 - “itfl 



3g 3{J 3w i 

3X^ _t °i3X^ 3X^ ~i 
(Eq. 24) 



The left hand side of the equation is zero because of Eq. 21, 



and thus 



3w f 

1 

3X„ 



,T 

*i 



3 K 



3M 



3X 



- to ? 



I 



i 3X 






(Eq. 25) 



d)T M <t> . 

*1 a ~1 
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The partial derivative of the generalized mass matrix with 
respect to the reciprocal of the area variable is defined as: 



3.M 



ax, 



:G = _ 



NE 

Z 

i=l 



p. AfL. 
ill 



i£ 



2I 3 Z 3 
Z 3 2I 3 



(Eo. 26) 



The partial derivative of the lumped mass matrix with respect 

to the reciprocal of the area variable is formulated following 

the same procedure as the generalized mass matrix. 

4 . Gradient of Stresses with Respect to the Joint 
Coordinate Variables 



Since stress in a member is 
K . u • 

»i y i 



0 . = 

~i 



A. 

l 



(Eq. 27) 



the gradient is calculated as 



35i 3 i n 3^ 

— = _L_ u . + K. — 

ax £ ax z -1 A. =1 ax £ 



(Eq. 28) 



In this case, the stiffness matrix is a function of the 
coordinate variables so the first term on the right hand 
side of Eq. 28 is not zero. 

5 . Gradient of Displacements with Respect to the Joint 
Coordinate Variables 



Consider the following equation: 
Ku = F 

The gradient of u with respect to X^ is: 
3 (u) = -K _1 K)u 



3X, 



3X £ 



(Eq. 29) 



(Eq. 30) 
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The partial derivative of the global stiffness matrix K now 
is defined as: 




NE 

I 

i=l 



A. E. 

l l 



% 

X 



3A . E . 

(-1) j (dx. ) [K. ] 
2 1 



Lf 



" (Eq. 31) 

where (-1)- 1 and j = 1 or 2 defines the sign of the gradient 
at that particular node. Note that the terms in the summa- 
tion are only evaluated for members connected to the joint 
defined by the variable X^. 

6 . Gradient of the Natural Frequencies with Respect to 
the Joint Coordinate Variables 

The value of the derivative of the eigenvalue w 2 
found from 



[K - u?M]<|>. = 0 

sj 1 ~ ~ 1 



(Eq. 32) 



is 
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(Eq. 33) 



The partial derivative of the generalized mass matrix with 
respect to the coordinate variables are defined as: 
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3X 
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- I 



(-1) 
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o . A . dX . 

l ii l 
6L i 



2I 3 I 3 
I 3 2I 3 



(Eq. 34) 



where (-l) 3 and j = 1 or 2 defines the sign of the gradient 
at the particular node. The partial derivative of the lumped 
mass matrix follows the same procedure. 
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D. APPROXIMATION CONCEPTS 



At this stage the analysis tools necessary for the design 
process under static and dynamic conditions may be regarded 
as available. It has been pointed out that the application 
of mathematical programming methods for structural design is 
the most widely used because of its great generality and its 
simple formulation. However, it is required that a large 
number of structural analyses and sensitivity analyses be 
performed. This has motivated the idea of formulating simple 
and explicit approximations for the most relevant response 
quantities . 

These approximations can only be expected to be of 
acceptable quality in some finite region of the design space 
surrounding the point about which the approximations were 
constructed. The total number of analyses required to find 
an optimum design using approximation concepts is significant- 
ly less than the number previously required. 

Some sources of simplification can be considered. First, 
the dimension of the design space can be reduced if a proper 
subspace can be identified. Second, linking of design 
variables which is imposed because of symmetry or practical 
considerations also reduces the dimension of the design 
space . 

The objective function for the design of trusses is a 
relatively simple explicit function of the design variables. 
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On the other hand nonlinear constraint functions are very 
complicated. From a computational standpoint, a small 
portion of the constraints play an active role in the opti- 
mization process; therefore, deletion of non-critical 
constraints avoids effort of evaluation of irrelevant 
constraints . 

A method to deal with complicated constraint which is 
very effective in reducing computational effort is to deduce 
simple and explicit expressions for the constraints. 
Linearization is directly and efficiently accomplished by 
Taylor series expansions. The order of the expansion 
selected is decided based on the degree of nonlinearity and 
the approximation required; and a trade-off must be made 
between the computational effort required for the highest 
order derivatives versus improvement of approximation. 

Application of a first order Taylor series expansion on 
stress constraints has been found to be sufficient. This 
was not the case for the natural frequency constraints, 
which may still be numerically unstable during the optimi- 
zation. This suggests further expansion up to a second 
order. However, the computational effort needed to do this 
usually exceeds the benefits. Therefore, a first order 
expansion with move limits is used. 

First order Taylor series expansion of a function, W, 
of the variable, X, about a point, Xo, is written as 
W(p)s W(x) = W(Xo) + (X-Xo) VW(Xo) 
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The philosophy underlying the use of linear approxima- 
tions is not to transform the problem into a sequence of 
linear programs, but to replace the constraints by simple 
and explicit approximation functions. In this study the 
Taylor series expansion was used for all constraints with 
respect to the reciprocal of the member sizing variables, A^ . 
The objective function, in reciprocal space, is now non- 
linear but still explicit and easily evaluated. 

E. OBJECTIVE FUNCTION 

The function whose least value is sought in the optimi- 
zation procedure is defined as the total weight of the truss 
which is given by 

NE 

W = Y p.A.L. (Eq. 35) 

JL . . x 1 1 
1=1 

Where NE is the total number of members, and is the 
material density, both are prescribed constants. A^ and 
L i are the area and length of the ith member respectively. 

F. CONSTRAINTS 

The restrictions to be satisfied in order for the design 
to be acceptable are formulated explicitly; behavioral and 
side constraints are defined accordingly. 

1 . Stress 

a . < a . < cu . (Eq. 36) 

ci — 1 — ti 
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where 0 c ^ and 0 ^ are the allowable compressive and tensile 
stresses respectively for member i. 

2 . Euler Buckling 



The Euler buckling stress is formulated as 



cA.E. 
— — 11 

°bi " L? 



(Eq. 37) 



where c is a prescribed constant. 

3 . Frequency 



io T , ^ w T , ^ u> 



K 



K 



(Eq. 38) 



where u> K and w are the lower and upper bounds respectively 
on the first natural frequency of the system. 

4. Limits on Areas 

(Eq. 39) 



A . 5 A. 1 A 

mm^ 1 max^ 



where, A . , and A are the minimum and maximum allowable 

mm . max . 

1 1 

cross-sectional areas of the ith member, and are taken to 
be the same for all members. When symmetry of design is to 
be preserved, linking of the variables is required. This 
is defined as 

Aj, = A i (Eq. 40) 

where K and i are symmetric members. 

5 . Displacement 



% “ ~ U K 



(Eq. 41) 



where u„ and u^ are the lower and upper bounds respectively 

K K 



on joint displacements u 



K' 
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6 . Limits on Coordinate Variables 



X K “ k ~ (Eq. 42) 

where again and X^ are the lower and upper bounds 
respectively on the Kth coordinate variable, and if symmetry 
must be preserved, linking of the variables is required. 

This is obtained by 

~ a l + kjl X K (Eq* 43) 

where a^ and b^ are constants and X^, is the coordinate 
variable . 



G. GENERAL FORMULATION 

The inequality constraints are of the form, 

j = 1 NIC 

The constraints are normalized as follows: 



Gj (X) < 0 



(Eq. 44) 



Stresses 



— - 1 - 0 



rr- -1-0 

a c 



(Eq. 45) 



Euler Buckling S- -1^0 (Eq. 46) 

<?b 

Note that because a and depend on the member area A^ , 
this constraint is treated as a nonlinear function. The 
Taylor series expansion is performed on o and the value 
of is continually updated. 
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Displacements 



— - 1-0 
6 - 



Frequency 



— - 1 - 0 

6 + 

3 - 1-0 

u 



The side constraints are 



A £ ' A i - ° 



X K - X K - ° 
X K - X K 1 0 



(Eq. 47) 



(Eq. 49) 



(Eq. 50) 
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III. OPTIMIZATION 



A. INTRODUCTION 

The main goal of structural engineering is to design 
structural systems that efficiently perform specified 
purposes. Selection of a specific algorithm is required 
and this algorithm must minimize the number of times the 
structure has to be analyzed and the amount of specific 
gradient information required. Finally, the algorithm 
should provide reasonable assurance that it will attain an 
optimum or near-optimum design. 

The next two sections are a brief explanation of the 
algorithm used for this work. 

B. GENERAL FORMULATION 

The general constrained minimization design problem is 
defined as 
minimize 

W(X) (Eq. 1) 

subject to 

Gj (X) - 0 j = 1, .m (Eq. 2) 

where, W(X) is the objective function. Functions G^ (X) are 
the set of inequality constraints. The vector of design 
variables X includes member sizing variables X^ and geo- 
metric design variables X^.'. 
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C. OPTIMUM GEOMETRY DESIGN 

The procedure used here was to treat the geometric design 
parameters as independent design variables. The member sizing 
parameters are handled as dependent variables which are 
determined as a sub-problem. 

Beginning with an initial geometric design vector X Q , the 
design proceeds iteratively using the following relationship: 



Xf 1 = X q + a* S q 

g g g g 



(Eg. 3) 



where q is the iteration number and S is the search direction 

g 

* 

to be determined. is the scalar parameter determining the 

distance of travel in the design space. 

For each proposed geometric vector, X^, the structure 
was optimized with respect to the member sizing variables, 

X , by the sub-optimization problem defined in Section D. 

Assume that for the initial geometry the structure has 
been optimized with respect to the cross-sectional areas, 
and that, from this subproblem there are £ active constraints 
of the form: 

Gj (X) = 0 k = 1, . . . £ (Eq. 4) 

where, Gj (X) is defined as active if its value is close to 
zero . 

Now, it is necessary to find the search direction, S , 

y 

so that by moving in this direction in the coordinate 
design space, the objective function is minimized. This 
direction is found by solving the following subproblem. 



32 



Minimize 



VF (X) . S (Eg. 5) 

Subject to: 

G r (X) . S < 0 k = 1,....£ (Eq. 6) 

S . § * 1 (Eq. 7) 

At this point, S, provides a search direction in the 

combined [X + X^] space which is projected onto the active 

constraints. Only the S g part of s is used. 

Once, S , is known it is substituted back into Eq. 3 

and a one-dimensional search on is performed to update 

the vector, X . The Optimum-Geometry design problem can be 

y 

summarized in the following algorithm. 

1. Given an initial coordinate design vector, X Q , and 

area design vector, A Q . Specify, D max and P max . <D max , 

maximum change in the coordinate at each iteration, and 

P , total number of iterations.) 

max 

2. Solve the fixed geometry problem and calculate the 
minimum weight (W*) for the current geometry. 

3. Determine the set of active constraints. 

4. Determine the search direction vector, S. 

•k 

5. Find the move parameter a, for a D change of some 

1 max 

coordinate variable, or such that some coordinate constraints 
becomes active. 

6. Solve the fixed geometry problem for coordinates. 

X = X + a*S and calculate the optimum weight, W (X) . 
g g 1 g 
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, * 

7. Find the parameter a 2 , which minimizes the weight 

"Jc 

using 2-point quadratic interpolation. If « 2 > cu do to 
step 9 . 

8. Solve the fixed geometry design problem for coordi- 

-q+1 -!Q * — Q 

nates = X^ + c* 2 and calculate the optimal weight. 

9. Check convergence; if satisfied, terminate; otherwise, 
set q = q+1 , update X and return to step 3 . 

D. FIXED GEOMETRY DESIGN 

As stated earlier, a sub-optimization problem has to be 
solved for a proposed geometry design vector, X. The 
structure is now optimized based on the cross-sectional 
area subspace. This is defined in general as: 

Minimize, 

W(X) (Eq. 8) 

Subject to, 

Gj (X) < 0 j = 1, Z (Eq. 9) 

The design proceeds iteratively. Given an initial vector 
of design variables X Q , find the X vector at the (q+1) ith 
iteration defined as: 

X q+1 = X q + X* S q (Eq. 10) 

m m m m 

* . . - 
where, a , is a scaler multiplier, and, S , is a vector move 
m m 

direction in the design space. 

Now the problem becomes one of finding the direction, 

S, and the move parameter, a*. Zoutendijk [Ref. 12] shows 
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that the direction may be found by solving the following 
problem. 



Maximize 



8 



(Eq. 11) 



Subject to 

VF (X) . S + 8 - 0 



(Eq. 12) 



G j (X m } • S + 0j 8 S 0 



j = 1 



f • • • 



NAC. 



(Eq. 13) 



S . S < 1 



(Eq. 14) 




If some of the constraints are violated, this algorithm 
is modified in order to find a feasible design [Ref. 13] . 

For the fixed geometry sub-problem, the reciprocals of 
the member sizes are used as design variables and approxima- 
tion techniques are employed [Ref. 5]. The optimization 
program CONMIN [Ref. 14] is used to solve this sub-problem. 
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IV. NUMERICAL EXAMPLES 



A. INTRODUCTION 

Design of planar trusses and space towers are presented 
here and the corresponding numerical results are summarized 
to demonstrate the purpose of this study. 

The examples begin with an 18-bar truss. In this case, 
it is shown how the optimum geometry is dependent on the 
constraints imposed. A single load condition was considered 
and the cross-sectional areas are linked. 

Next, a 47-bar planar tower is designed to support a set 
of load conditions given in Table VI. The design is subject 
to constraints on the member stresses, Euler buckling, dis- 
placement and first natural frequency. Linking was imposed 
for symmetry in both cross-sectional area and coordinate 
variables . 

Finally, a 25-bar space tower was designed. Constraints 
on stresses, Euler buckling, displacement, and frequencies 
were imposed. The truss was required to support two 
different load conditions. Symmetry for the member areas 
and coordinates was established. 

B. CASE 1: 18-BAR PLANAR TRUSS 

A cantilever truss, as shown in Fig. A.l, has been used 
previously as an example for the design of trusses of a 
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specified geometry [Ref. 7]. This structure was designed 
for optimum geometry subject to a single set of load condi- 
tions given in Table I. The allowable stresses are speci- 
fied as 

-20,000 £ ck $ 20,000 psi 

Young's modulus is taken as 10^ psi and the material density 
P = 0.1 lb./cu. in. The allowable stress at which Euler 
buckling occurs is 



The independent coordinate variables were taken as X3, 

Y3, X5, Y5, X7, U7 , X9 , Y9 . The member areas were linked 
in the following groups: A1=A4=A8=A12=A16 ; A2=A6=A10=A14=A18 ; 

A3=A7=A11=A15 ; A5=A9=A13=A17 . There are a total of eight 
independent coordinate variables and four independent area 
variables . 

1 . Case la 

The structure was designed subject to stress constraints 
only. The resulting geometry is shown in Fig. la, and the 
design information is given in Table II. Weight versus 
iteration history is plotted in Fig. A. 2. The number of 
analyses for this design is 59, and the execution time 2.53 
seconds on an IBM 3033 computer. 
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2. Case lb 



Stress and Euler buckling constraints were imposed 
for this case. Design information is given in Table III. 
Figure lb shows the final geometry layout and the weight 
versus iteration history is shown in Fig. A. 3. The number 
of analyses for this design is 78 and the execution time 
was 3.94 seconds. 

3. Case 1c 

This design is based on stress, Euler buckling and 
displacement constraints. The latter was applied at node 
number 1 in the Y-direction. The results are shown in 
Table IV. Figure lc presents the final layout and the 
weight versus iteration is plotted in Fig. A. 4. The number 
of analyses is 91 and the execution time was 3.94 seconds. 

4 . Case Id 

This final case includes all the constraints mentioned 
before plus a constraint that the first natural frequency 
be greater than or equal to 3 Hz. A non-structural mass 
of fixed value W = 1,000 lbs was placed at node 1. Steady 
convergence is achieved and results are summarized in Table V. 
The final geometry and weight iteration history are shown in 
Figures Id and A. 5, respectively. The number of analyses is 
96 and the execution time was 4.73 seconds. 

When the mass is removed from the structure, the 
design fails to converge even when move limits are imposed. 
First, it is known by definition that the frequency is 
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proportional to the stiffness and is inversely proportional 
to the mass. Second, both the K and M matrices are functions 
of area and coordinate variables. Therefore, as K increases, 
M also increases, and the frequency is kept closed to the 
initial range. 



C. CASE 2: 47-BAR PLANAR TOWER 



The initial layout of the tower is shown in Fig. A. 6. 
Stress constraints were imposed as well as constraints on 
Euler buckling, displacement, and first natural frequency: 



and 



-15,000 < a i < 20,000 psi 




10. IttEA. 

i 




psi 



i=l ,47 



X * 12 hz. 



respectively. The members are assumed tubular with D/t=10. 

7 

Young's modulus of 3.x 10 psi, and material density 
£ = .3 lbs/cu. in. Minimum allowable area of 10 ^ in. was 
imposed. Symmetry about the y-axis was desired during 
optimization so linking of variables is necessary. Joints 
15, 16, 17, 22 are kept fixed, and joints 1 and 2 allowed 
to move in the x-axis direction (Y=0) . This gave a total 
of 27 area design variables and 17 coordinate variables. 
Nonstructural masses were attached at nodes 19 and 20, of 
W=500 lbs each. The results are tabulated in Tables VII 
through XII. The final geometry for the case where all 
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constraints were imposed is shown in Fig. A. 7 and the 
iteration history is shown in Figs. A. 8 through A. 10. When 
all constraints are imposed, the design required 220 analysis 
and the execution time was 160.7 seconds. 



D. CASE 3: 2 5 -BAR SPACE TOWER 



The 25-bar space tower shown in Fig. A. 11 was designed 
to support two independent load conditions given in Table 
XIII. The allowable stresses were specified as 
-40,000 < a i £ 40,000 psi i=l,25 

Young's modulus was selected as 10 7 psi and the material 
density p = .1 lb./cu. in. The members are assumed to be 
turbular with a nominal diameter to thickness ratio of 
D/t = 100., so that the stress at which Euler buckling occurs 



is 



bi 



10.1 7T EA. 

l 

8L? 



psi , 



Symmetry with respect to both x-z plane and the y-z plane 
was imposed, so linking of variables were made as follows: 
the member areas were grouped in the following sequence: 

Al, A2=A3=A4=A5, A6=A7=A8=A9 , A10=A11, A12=A13, A14=A15=A16=A17 , 
A13=A19=A20=A21, and A22=A23=A24=A25 . For the coordinates 
X4 , Y4, Z4, X8 , Y8 were considered as variables. The joints 
1 and 2 were held fixed and joints 7 through 10 were required 
to lie in the x-z plane. Nonstructural masses of W=500 lbs. 
were attached at nodes 1 and 2, respectively. The first 
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natural frequency was limited to a value \ ^ 16 Hz. 
Results of this example are shown in Table XIV and the 
iteration history is shown in Fig. A. 12. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The final layout dependency on constraints has been 
presented for the elastic design of trusses for optimal 
geometry. The truss may be planar or three-dimensional and 
may be indeterminate. Multiple load conditions were con- 
sidered. The design procedure was separated as: analysis, 

design for fixed geometry, design for optimum geometry. 

The displacement method for static analysis and the 
subspace iteration method for dynamic analysis were applied. 

The sequential optimization based on two design sub- 
spaces present substantial advantages in the reduced 
number of analyses and allow the designer to keep control 
of the optimization process. 

Several examples were considered. In every case, signifi 
cant weight reduction was efficiently achieved. Also, the 
geometries obtained appear quite acceptable from an 
aesthetic as well as structural point of view. 

The graphs of Weight vs. Iteration number show that an 
acceptable design can be achieved in few iterations. 
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B. RECOMMENDATIONS 



The following recommendations may be of theoretical 
and practical value. 

1. Weight due to the structure itself and other design 
dependent external forces can easily be taken into account 
and should be considered in future studies. 

2. Application should be made to reasonable sized 
structures such as offshore towers and long span roof 
trusses . 

3. The principles and the procedure described herein 
can also be used for optimal design of frames as well as 
other structures. 
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APPENDIX A 








Frequency. 



Figure A.1. 18-Ear Planar Truss. Stress, Euler Buckling, 
Displacement, and Frequency Constraint. 
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Iteration number 



18-BRR PLRNRR TRUSS 
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Figure A. 3. 18-Bar Planar Truss. Height vs. Iteration Num- 
ber. Stress, Euler Buckling. 
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Figure A. 4. 18-Bar Planar Truss. Weight vs. Iteration Num- 
ber. Stress, Euler Buckling, Displacement. 
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Figure A. 5. 10-Bar Planar Truss. Weight vs. Iteration 
ber. Stress, Euler Buckling, Displacement 
Freguency. 
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Figure A. 6. 47-Bar Planar Tower. 
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Figure A. 7. 47-Ear Planar lower. Stress, Euler Buckling, 
Displacement, Frequency. 
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Figure A. 8. 47-Ear Planar Tower. Height vs. Iteration Num- 
ber. Stress, Euler Buckling. 
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Figure A. 9. 47-Ear Planar lower. Height vs. Iteration Num- 
ber. Stress, Euler Buckling, Displacement 
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Figure A. 10. 47-Bar Planar Toner. Height vs. Iteration Num- 
ber. Stress, Euler Buckling, Displacement, 
Freguency . 
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Figure A. 11. 25-Bar Space Tower. 
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25-BRR SPRCE TOWER 
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Figure A. 12. 25-Ear Space lower. Height vs. Iteration Num- 
ber. Stress, Euler Buckling, Displacement, 
Freguency . 



TA3LE I 



18-Bar Elanar Truss. Loads and Constants. 



Joint 

1 

2 

4 

6 

8 



LCAC CONDITION 
Fx 

0 . 

0 . 

0. 

0. 

0. 



(lbs.) 

Fy 

- 20 , 000 . 

- 20 , 000 . 

- 20 , 000 . 

- 20 , 000 . 

- 20 , 000 . 



TAELE CF CONSTANTS 



Young's Modulus 


E= 


0. 1EC8 


psi. 


(.69E08 


KN/m ) 




Allowable Stress 


O' = 


0.2E05 


Fsi. 


(. 1 38E06 


KN/m 


) 


Density 


r- 


0. 1 lbs 


/cu .in 


(. 276E04 


Kg/m 


) 


Buckling Coefficient 


C= 


. 4 EOS 


psi . 


(. 276E0S 


KN/m 


) 
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TABLE II 



* Areas 



18-Bar Planar Truss. Design Infor m ation . Stress. 



AREA* (sg.in.) 



flember 


Initial 


Final 


A 1 = A4 = A8 = A12=A1 6 


10. 0 


11.05 


A2=A6 = A10 = A14 = A 18 


15. 0 


15.07 


A3=A7=A1 1=A 15 


5.0 


4.54 


A4=A9=A13=A17 


7. 07 


5.33 



COORDINATES (in.) 



Joint 



Initial 



Final 



are 





X 


y 


X 


y 


3 


1000 


0.0 


99 1. 17 


19.686 


5 


750 


o 

• 

o 


745.88 


15. 12 


7 


500 


0,0 


494.61 


34. 538 


9 


250 


0.0 


249. 54 


23.609 



Optimum Height for Initial Geometry = 4,780.5 (Its.) 

Optimum Height for Final Geometry = 4,524.7 (lbs.) 

the optimum values for the initial and final geometry. 
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TAELE III 



18-3ar Planar Iruss. Design Information. Stress, Euler 

Buckling. 



A a E A * (sq.in.) 



Member 


Initial 


Final 


A 1=A4=A8=A12=A16 


10. 00 


11 . 


34 


A2=A6=A10=A14=A18 


21.65 


19, 


28 


A3=A7=A1 1=A15 


12. 5 


1C. 


97 


A 4= A9 = A 1 3 = A 1 7 


7.071 

CCOEDI NATES (in.) 


C. 


306 


Joint 


Initial 


Final 




x y 


X 


y 


3 


1000 0.0 


994.57 


162. 31 


5 


750 0.0 


747.36 


102.92 


7 


500 0.0 


482.90 


32. 962 


9 


250 0.0 


221.71 


17. 105 


Optimum Weight 


for Initial Geometry 


= 6,430.7 


(lbs. ) 


Optimum Weight 


for Final Geometry 


= 5,713.0 


(lbs. ) 



Areas are the optimum values for the initial and final geometry. 



TABLE IV 



13-Bar Planar Truss. Design Information. Stress, Euler 

Buckling, Displacement. 



AREA* (sg.in.) 



Member 


Initial 


Final 


A4 =A8 = A 1 2 = A 1 6 


1«. 18 


16.27 


A6=A1 0=A 1U=A 18 


21.66 


2 0.06 


A7 = A 1 1 = A 15 


12. 5 


11.18 


A9=A13=A 17 


10.32 


7.863 







COORDINATES (in.) 






Joint 




Initial 


Final 






x I 


X 


y 


3 




1000 0.0 


962.19 


109. 93 


5 




750 0.0 


703.95 


55.309 


7 




500 0.0 


452.7 


27. 307 


9 




250 0.0 


208.48 


- 2.884 


Optimum 


Height 


for Initial Geometry 


= 7,413. 1 


(lbs. ) 


Optimum 


Heig ht 


for Final Geometry 


= 6,466.9 


(lbs. ) 



* Areas are the optimum values for the initial and final gcemetry. 
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TABLE V 



18-Bar Planar Truss. Design Information. Stress, Euler 
Euckling, Displacement, Frequency. 



area* (sg.in.) 



Bember 


Initial 


Final 


A 1= A4=A8= A 1 2 = A 1 6 


31. 68 


29. 


66 


A2=A6=A10=A14=A 18 


40.58 


35. 


29 


A3= A7=A1 1=A 15 


12. 5 


13. 


03 


A 4 = A9=A 1 3=A 1 7 


21. 10 

COORDINATES (in.) 


16. 


09 


Joint 


Initial 


Final 




x y 


X 


y 


3 


1000 0.0 


972.52 


88. 0 


c 


750 0.0 


717.08 


27. 8 4 


7 


500 0.0 


468.75 


- 3.91 


9 


250 0.0 


200. 98 


-19.29 


Cptimun Weight 


for Initial Geometry 


= 13,686 . C 


(lbs. ) 


Optimum Weight 


for Final Geometry 


= 11 , 528. C 


(lbs. ) 



* Areas are the optimum values for the initial and final geometry. 
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TABLE VI 



47-Bar Planar Tower. Load Conditions. 



LOAE CONDITION 1 (lbs.) 



JOINT 




Fx 


Fy 


17 




6, COO. 


-14,000. 


22 




0. 


0. 


LOAD 


CONDITION 


2 


(Its .) 


Joint 




Fx 


Fy 


17 




0. 


0. 


£ t 


6,000. 


1 

% 

O 

o 

o 

• 


LOAD 


CONDITION 


3 


(its.) 


Joint 




Fx 


Fy 


17 




6, COO. 


-14,000. 


22 




6 ,0 00 . 


-14,000. 
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TAELS VII 



4 7- Bar Planar lexer. Area. Stress, Euler Euckling. 



Heater 


Initial 


Final 




Area * (sg.in.) 


Area 


3 


3.76U 


2.727 


4 


3.315 


2.468 


5 


0.787 


0.727 


7 


0.864 


0.213 


3 


0.856 


0. 938 


10 


1.754 


1.076 


12 


2.087 


1.691 


14 


1.188 


0 . 6 55 


15 


1 .52S 


1.058 


18 


2.087 


1.412 


20 


0.648 


0.263 


22 


0.843 


0.811 


24 


1 .700 


1.060 


26 


1 .700 


1.052 


27 


1.354 


0.820 


28 


0.847 


0.302 


30 


3.609 


2.766 


31 


1.435 


0.657 


33 


0.638 


0.207 


35 


2.842 


2. 657 


36 


0.676 


0.266 


38 


1.596 


1.408 


40 


3.686 


3.429 


41 


1.526 


0.991 


43 


0.677 


0.170 


45 


4.486 


3.650 


46 


1.532 


1.005 



* Cptimum areas for initial and final geometry. 
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47-Bar Planar lower. Coordinates. Stress, Euler Bucling 
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47-Bar Planar Tower. Coordinates. Stress, Euler Buckling 
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TABLE X 



47-3ar Plana; 



flemfcer 

3 

4 

5 

7 

8 

10 

12 

14 

15 
18 
20 
22 
24 
26 

27 

28 

30 

31 
33 

35 

36 
38 

40 

41 
43 

45 

46 



Tcwer. Areas. Stress, Euler Buckling, 



Displacement . 




Initial 


Final 


Area * (sq.in.) 


Area * (sq 


15.28 


4 .478 


16. 55 


3.671 


3.066 


1.037 


2.691 


0.561 


2.435 


1 . 156 


3.706 


1 . 170 


5.263 


2.318 


4.048 


0.669 


5.077 


1.385 


4.780 


1 . 169 


1 .650 


0.562 


2.231 


1.295 


4.665 


1 .526 


4.816 


1.511 


3.693 


1.036 


2.923 


0.566 


17.91 


5.231 


7.091 


0.949 


2.566 


0.558 


17.59 


6. 175 


6.640 


1.217 


3.385 


1.076 


28.1 1 


7.730 


6.792 


1.S70 


2.701 


0.555 


37. 1 1 


9.878 


7.363 


2.153 



* Optimum ares for initial and final geometry. 
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TABLE XI 



47-Bar Elanar Tower. Areas. Stress, Euler Buckling, 
Displacement, and Frequency. 



Heater 


Initial 
Area * 


3 


16.59 


4 


15.37 


5 


1.321 


7 


1 .871 


8 


1.356 


10 


5.152 


12 


8.593 


14 


5.988 


15 


8.468 


18 


8.593 


20 


1.416 


22 


1 .836 


24 


8.470 


26 


8.470 


27 


5.999 


28 


1 .949 


30 


14.54 


31 


7.096 


33 


1.709 


35 


9 .710 


36 


1.414 


38 


1.385 


40 


10. 43 


41 


1 .484 


43 


1.465 


45 


11.17 


46 


1.482 



Final 

(sq.ir.) Area * (sq.in.) 

4.625 
4.605 
0.S38 
0.348 
1.081 
1.089 
3.025 
1 .088 
1.691 
1 .656 
0.635 
1.724 
2.480 
2.186 
1.342 
0.452 
4.833 
0.855 
0.307 
4.579 
0.368 
1.439 
4.396 
1.126 
0.270 
4.365 
1.664 



* Optimum areas for initial and final geometry. 
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47-Bar Planar Tower. Coordinates. Stress, Euler Euckling 
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Optimum Weight for Final Geometry 



TABLE XIII 
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TABLE XIV 



25-Ear Space lower. Design Information. 







COORDINATES (in 


‘> 






Joint 




Initial 




Strs-fcck . 


Displaceaen t 


Frequency 




z 


37.5 




21.487 


33.512 


64.913 


4 


y 


37.5 




46.271 


602. 83 


47.271 




z 


100. 




ICO. 27 


1 14.26 


128.77 


8 


X 


100. 




22.146 


47.83 


107.56 


y 


100. 




96.353 


140.53 


110.077 








AREAS 


* (sq.in.) 






EENBER 




STRESS- 

INITIAL 


BUCKL 

FINAL 


DISPLACEMENT FREQUENCT 

INITIAL FINAL INITIAL FINAL 


1 




0. 9E-3 


0.013 


0.668 


0. 029 0.54b 


0.Q7O 


2 




0. 782 


0.414 


2.032 


0.450 2.533 


0.73b 


6 




0.754 


0.842 


3. 100 


1.219 2.635 


1.209 


10 




0. IE-3 


0.033 


0. 533 


0.015 0.547 


0.054 


12 




0. 130 


0.101 


0. 549 


0.118 C.496 


0.206 


14 




0.558 


0.121 


0.668 


0.084 0.533 


0.350 


18 




0.982 


0.739 


1.592 


0.751 1.449 


0.720 


22 




0.801 


0.554 


2.653 


0.668 3.363 


1. 130 





Stress-Buckl. 
Initial Pinal 


Displaceaent 
Initial fianal 


Frequency 
Initial Final 


Weiqht (lbs.) 


229.52 128.55 


565.2 


169.9 


587.73 


261.52 


Iteration 


12 




44 




9 


Ti»€ (sec.) 


5.38 




37. 14 




21.09 



* Optimuo areas for initial and final geoietry. 
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F2533 Felix 

c.l Shape optimization 

of trusses subject to 
strength, displacement, 
and frequency con- 
straints. 



